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ABSTRACT 

The recent discovery of a high optical depth r to Thomson scattering from the Wilkinson Microwave Anisotropy 
Probe (WMAP) data implies that significant reionization took place at redshifts z > 6. This discovery has important 
implications for the sources of reionization, and allows, for the first time, constraints to be placed on physical 
reionization scenarios out to redshift z ^ 20. Using a new suite of semi-analytic reionization models, we show 
that the high value of r requires a surprisingly high efficiency e of the first generation of UV sources for injecting 
ionizing photons into the intergalactic medium. We find that no simple reionization model can be consistent with 
the combination of the WMAP result with data from the z ^ 6.5 universe. Satisfying both constraints requires either 
of the following: (i) H2 molecules form efficiently at z ^ 20, survive feedback processes, and allow UV sources 
in halos with virial temperatures Tvir < lO'^K to contribute substantially to reionization, or (ii) the efficiency e in 
halos with T^ii- > 10*K decreased by a factor of > 30 between z 20 and z ^ 6. We discuss the relevant physical 
issues to produce either scenario, and argue that both options are viable, and allowed by current data. In detailed 
models of the reionization history, we find that the evolution of the ionized fractions in the two scenarios have 
distinctive features that Planck can distinguish at <; 3a significance. At the high WMAP value for r, Planck will 
also be able to provide tight statistical constraints on reionization model parameters, and elucidate much of the 
physics at the end of the Dark Ages. The sources responsible for the high optical depth discovered by WMAP 
should be directly detectable out to z ~ 15 by the James Webb Space Telescope. 



1. INTRODUCTION 

How and when the intergalactic medium (IGM) was reion- 
ized is one of the long outstanding questions in astrophysical 
cosmology, holding many clues about the onset of structure for- 
mation in cold dark matter (CDM) cosmologies, and the nature 
of the first generation of light sources. The past year has seen 
an explosion of progress in probing the reionization history, 
culminating in recent results from the WMAP satellite (Ben- 
nett et al. 2003). Spectroscopic observations of high redshift 
quasars had previously found no strong HI absorption (a so- 
called "Gunn-Peterson ", GP trough) in sources at z ^ 6, show- 
ing that the IGM is highly ionized at these redshifts (Fan et al. 
2000). However, the discoveries of a handful of bright quasars 
in the Sloan Digital Sky Survey (SDSS) at redshifts z > 6 have 
revealed full GP troughs, i.e. spectra consistent with no flux 
at high S/N over a substantial stretch of wavelength shortward 
of (1 +z)A„ = 8850A (Becker et al. 2001; Fan et al. 2003). 
The lack of any detectable flux translates to strong lower limits 
xh ^ 0.01 on the mean mass-weighted neutral fraction of the 
IGM at z - 6 (Cen & McDonald 2002; Fan et al. 2002; Lidz et 
al. 2002; Pentericci et al. 2002). Comparisons with numerical 
simulations of cosmological reionization (Gnedin 2001; Cen & 
McDonald 2002; Fan et al. 2002), together with the rapid rise 
towards high redshifts of the neutral fractions inferred from a 
sample of high redshift quasars from 5 < z Ss 6 (Songaila & 
Cowie 2002), suggest that the IGM becomes neutral at a rela- 
tively narrow redshift interval, Az ~ 1 , beyond the redshifts of 
the most distant sources. 

An alternative way of probing the reionization history is to 
study the temperature and polarization anisotropics of the cos- 
mic microwave background (CMB). Both anisotropy patterns 
are affected by Thomson scattering of CMB photons off the 
free electrons in a reionized IGM (Zaldarriaga 1997). Physi- 
cally, the CMB and the GP troughs probe two different stages 



of reionization: the CMB is sensitive to the initial phase when 
Xh first decreases below unity, and free electrons appear; say, 
at redshift Ze- On the other hand, the (hydrogen) GP trough is 
sensitive to the late stages, when the remaining neutral hydro- 
gen atoms are cleared away, at redshift zh (Kaplinghat et al. 
2003). The recent measurement by the WMAP (Bennett et al. 
2003) satelhte of the optical depth, r = 0.17 ±0.04 to Thomson 
scattering (Kogut et al. 2003) implies a reionization redshift 
of Ze = 17 ±5, under the assumption that reionization occurs 
abruptly. This is to be compared to zh ^6-1, obtained from 
extrapolations of the available high redshift quasar spectra at 
z ~ 6. If the redshift zh coincided with Ze, it would imply opti- 
cal depths of only r ~ 0.4-0.5. While the formal discrepancy 
between this opacity and the WMAP value is only at the ^ 3<j 
level, it is timely to ask the following questions: (1) what are 
the physical requirements on the ionizing sources in a ACDM 
cosmology to produce an optical depth ofr ^ 0.17, and (2) are 
these consistent with other data at z ^^,6? 

In this paper, we address both of these questions. In the 
ACDM cosmology favored by WMAP, ~ 3 sigma peaks of the 
primordial density field on mass scales of M ~ 10^ M0 col- 
lapse at redshifts of z ~ 20. The onset of reionization at similar 
redshifts is therefore dependent on efficient formation of ioniz- 
ing sources in nonlinear objects down to these low mass-scales. 
Early reionization can be accommodated in models of reioniza- 
tion in ACDM cosmologies (Shapiro, Giroux & Babul 1994; 
Tegmark et al. 1994; Haiman & Loeb 1997, 1998; Valageas & 
Silk 1999; Haiman, Abel & Rees 2000), and is indeed a natural 
consequence in models including a sufficient number of metal- 
free stars (Wyithe & Loeb 2003 and Cen 2003). If the optical 
depth is as high as the central value determined by WMAP, then 
future CMB polarization measurements can provide significant 
constraint on such models (Kaplinghat et al. 2003, see also 
Venkatesan 2000, 2002). 

This paper is organized as follows: in § 2, we discuss the 
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relevant physical effects that can lead to reionization histories 
consistent with both the WMAP data and the z ~ 6 quasar spec- 
tra . In § 3, we then describe a semi-analytical model that takes 
these effects into account, and in § 4, we derive and discuss 
the ionization histories in a suite of models. In § 5, we quan- 
tify constraints on the models that will be available from future 
polarization measurements by WMAP and by the Planck satel- 
lite.' Motivated by the required presence of ionizing sources 
at z ~ 20, in § 6 we consider the highest redshift at which the 
James Webb Space Telescope can delect these sources. Finally, 
we discuss our results and offer our conclusions in § 7. In 
this paper we adopt the background cosmological parameters 
as measured by the WMAP experiment (Spergel et al. 2003, Ta- 
bles 1 and 2), n„ = 0.29, f^A =0.71, nh = 0.047, h = 0.72 and 
an initial matter power spectrum P(k) oc A;" with n = 0.99 and 
normaUzation ug = 0.9; with the exception of § 4.6, where we 
consider the running index model favored by the combination 
of WMAP with other data (Spergel et al. 2003, Table 10). 

2. SUMMARY OF RELEVANT PHYSICS AT HIGH REDSHIFTS 

In the context of ACDM cosmologies, it is natural to iden- 
tify the first dark matter condensations ("halos") as the hosts 
of the first sources of light. A simplified picture for reion- 
ization is as follows: the gas in dark matter halos cools and 
turns into ionizing sources, which drive expanding ionized re- 
gions into the IGM. The volume filling factor of ionized regions 
then tracks the formation of dark halos. Eventually, the ionized 
regions percolate, and the remaining neutral hydrogen in the 
IGM is quickly cleared away as the ionizing background builds 
up. From the point of view of the large-angle CMB polariza- 
tion anisotropics, what matters most is the free electron fraction 
Xe = rie/ritot- Accordingly, we now briefly describe the various 
effects that should determine the evolution of Xe(z) at high red- 
shifts. For in-depth discussion, we refer the reader to extended 
reviews (Barkana & Loeb 2001 ; Loeb & Barkana 2001), and to 
reviews focusing on the roles of H2 molecules in reionization 
(Abel & Haiman 2000), on the effect of reionization on CMB 
anisotropics (Haiman & Knox 1999), and on progress in the last 
two years in studying reionization (Haiman 2003). 

In the discussions that follow, it will be useful to distinguish 
dark matter halos in three different ranges of virial tempera- 
tures, as follows: 

lOOK ^ Tvir ^ lO^K (Type II) 

lO^K ^ Jvir ^ 2 X lO^K (Type la) 
Tvir ;^2x lO^K (Type lb) 

We will hereafter refer to these three different types of halos 
as Type II, Type la, and Type lb halos. The motivation in dis- 
tinguishing Type II and Type I halos is based on H2 molecular 
vs atomic H cooling (§ 2.1), whereas types "a" and "b" reflect 
the abiUty of halos to allow infaU and cooUng of photoionized 
gas (§ 2.3). While the actual values of the temperatures, and 
the sharpness of the transition from one category to the next is 
uncertain, as we argue below, each type of halo likely plays a 
different role in the reionization history. In short. Type II ha- 
los can host ionizing sources only in the neutral regions of the 
IGM, and only if H2 molecules are present; Type la halos can 
form new ionizing sources only in the neutral IGM regions, but 

' CMB constraints are discussed in more detail in a companion paper (paper 
n, Holder et al. 2003), focusing on information beyond measuring r, and on 
the possible bias in the values of t measured without a-priori knowledge of the 
reionization history. 



irrespective of the H2 abundance, and Type lb halos can form 
ionizing sources regardless of the H2 abundance, and whether 
they are in the ionized or neutral phase of the IGM. In the fol- 
lowing sections, we provide a summary of the various impor- 
tant physical effects. Much of this summary relies on existing 
hterature; our primary difference from existing ideas is in the 
description of the transition from metal-free to normal stellar 
population (§ 2.3). 

2.1. ¥l2 Molecule Formation 

In ACDM cosmologies, structure formation is bottom-up: 
the earliest nonlinear dark matter halos form at low masses. 
Gas contracts together with the dark matter only in dark halos 
above the cosmological Jeans mass, Mj « lO''. However, this 
gas can only cool and contract to high densities in somewhat 
more massive halos, with M ;> Mm = 1O^M0[(1 +z)/ll]-^/2 
(i.e. Type II halos), and only provided that there is a sufficient 
abundance of H2 molecules, with a relative number fraction at 
least nm/nH ~ 10"^ (Haiman, Thoul & Loeb 1996; Tegmark 
et al. 1997). Because the typical collapse redshift of halos is 
a strong function of their size, the abundance of H2 molecules 
is potentially the most important parameter in determining the 
onset of reionization. For example, 2a halos with virial tem- 
peratures of lOOK appear at z = 16, while 2a halos with virial 
temperatures of 10"*K (Type la halos, in which gas cooling is 
enabled by atomic hydrogen lines) appear only at z = 9. Table 1 
summarizes the collapse redshifts and masses of dark halos at 
these two virial temperatures. As a result, the presence or ab- 
sence of H2 in Type II halos makes a factor of r^2 difference in 
the redshift for the onset of structure formation, and thus poten- 
tially effects the reionization redshift and the electron scattering 
opacity r by a similar factor In the absence of any feedback 
processes, the gas collecting in Type II halos is expected to be 
able to form the requisite amount of H2 (Haiman, Thoul & Loeb 
1996; Tegmark et al 1997). However, both internal and exter- 
nal feedback processes can alter the typical H2 abundance (see 
discussion in § 2.5 below). 

2.2. The Nature of Ionizing Sources 

A significant uncertainty is the nature of the ionizing sources 
turning on inside halos collapsing at the highest redshifts. Three 
dimensional simulations using adaptive mesh refinement (AMR; 
Abel, Bryan & Norman 2000, 2002) and smooth particle hydro- 
dynamics (SPH; Bromm, Coppi & Larson 1999, 2002) tech- 
niques have followed the contraction of gas in Type II halos at 
high redshifts to high densities. These works have shown con- 
vergence toward a temperature/density regime of T ~ 200 K, 
n ~ lO'^cm"^, dictated by the critical density at which the ex- 
cited states of H2 reach equilibrium population levels. These 
results have suggested that the first Type II halos can only form 
unusually massive stars, with masses of at least ~ 200 M0, and 
only from a small fraction ( ^ 0.01) of the available gas. Such 
massive stars are effective producers of ionizing radiation, mak- 
ing reionization easier to achieve. In addition, it is expected 
that the first massive stars, forming out of metal-free gas, have 
unusually hard spectra (Tumlinson & ShuU 2000; Brormn, Ku- 
dritzki & Loeb 2001; Schaerer 2002), which is important for 
possibly ionizing helium, in addition to hydrogen. An alter- 
native possibility is that a similar fraction of the gas in the first 
Type II halos forms massive black holes (Haiman & Loeb 1997; 
Bromm & Loeb 2002). These early black holes can then ac- 
crete gas, acting as "miniquasars", and producing a hard spec- 
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tmm extending to the soft X-rays. These soft X-rays could be 
important in catalyzing H2 formation (see discussion in § 2.5 
below). 

2.3. Efficiency Parameters and the Transition from Metal Free 

to "Normal " Stars 

Another fundamental question of interest is the efficiency at 
which the ionizing source(s) associated with each halo eject 
ionizing photons into the IGM. This can be parameterized by 
the product e* = A'-,-/* /esc, where /* = M*/(nb^haio/f^m) is the 
fraction of baryons in the halo that turns into stars; A^^ is the 
mean number of ionizing photons produced by an atom cycled 
through stars, averaged over the initial mass function (IMF) of 
the stars; and /esc is the fraction of these ionizing photons that 
escapes into the IGM. It is difficult to estimate these quanti- 
ties at high-redshifts from first principles, but the numbers and 
discussion below can serve as a useful guide. 

Although the majority of the baryonic mass in the local uni- 
verse has been turned into stars (Fukugita, Hogan & Peebles 
1998), the global star formation efficiency at high redshifts was 
likely lower. To explain a universal carbon enrichment of the 
IGM to a level of 10"^- 10"^ Z©, the required efficiency, aver- 
aged over all halos at z ^ 4, is /* = 2-20% (Haiman & Loeb 
1997). However, the numerical simulations mentioned in § 2.2 
above suggest that the fraction of gas turned into massive stars 
in Type II halos is /* ^ 1%. 

The escape fraction of ionizing radiation in local starburst 
galaxies is of order ~ 10%. The higher characteristic densities 
at higher redshifts could decrease this value (Dove et al. 2000; 
Wood & Loeb 2000), although there are empirical indications 
that the escape fraction in z ^ 3 galaxies may instead be higher 
(Steidel, Pettini & Adelberger 2001). Radiation can also more 
readily ionize the local gas, and escape from the small Type II 
halos which have relatively low total hydrogen column densi- 
ties ( <, 10^^ cm-2). 

The ionizing photon yield per proton for a normal Salpeter 
IMF is A'^ K. 4000. However, if the IMF consists exclusively 
of massive M ^ 200 M0 metal-free stars, then A^^ can be up 
to a factor of ^ 20 higher (Bromm, Kudritzki & Loeb 2001; 
Schaerer 2002). The transition from metal-free to a "normal" 
stellar population is thought to occur at a critical metalUcity of 
Zcr ^ 5 X 1O~"*Z0 (Bromm et al. 2001). It is natural to associate 
this transition with that of the assembly of halos with virial tem- 
peratures of > lO'^K (Type la halos) . Type II halos are fragile, 
and Ukely blow away their gas and "shut themselves off" af- 
ter a single episode of (metal-free) star-formation. They are 
therefore unlikely to allow continued formation of stars with 
metallicities above Zoit. Subsequent star-formation will then 
occur only when the deeper potential wells of Type la halos are 
assembled and cool their gas via atomic hydrogen Unes. The 
material that collects in these halos will then have already gone 
through a Type II halo phase and contain traces of metals. 

As we argue in § 2.5 below, there exists an alternative, equally 
plausible scenario. Most Type II halos may not have formed 
any stars, due to global H2 photodissociation by an early cos- 
mic soft-UV background. In this case, the first generation of 
metal-free stars must appear in Type la halos. Halos above this 
threshold can eject most of their self-produced metals into the 
IGM, but, in difference from Type II halos, can retain most of 
their gas (MacLow & Ferrara 1999), and can have significant 
episodes of metal-free starformation. These halos will also 
start the process of reionization by driving expanding ioniza- 



tion fronts into the IGM. The metals that are ejected from Type 
la halos will reside in these photoionized regions of the IGM. 
As discussed in § ref:feedback below, photoionization heating 
in these regions suppresses gas infall and cooling, causing a 
pause in the formation of new structures, until larger dark mat- 
ter halos, with virial temperatures of Jvir 2 x lO^K (Type lb 
halos) are assembled. The material that collects in Type lb ha- 
los will then have already gone through a previous phase of 
metal-enrichment by Type la halos, and it is unlikely that Type 
lb halos can form significant numbers of metal-free stars. 

The effect of a metal-free to a normal stellar population on 
the reionization history was studied recently by Wyithe & Loeb 
(2003) and Cen (2003). Both of these works assumed that 
rather than being tied either to a Type II— +Ia or Type la^Ib 
halo transition, the switch from metal-free to normal stars oc- 
curs abruptly at a critical redshift Zcr- This would be justified 
in the picture of a uniformly increasing metal-enrichment of 
the IGM, and would lead to qualitatively different reionization 
histories from the ones we obtain below. 

2.4. Density Fluctuations ( dumpiness) of the IGM 

As the ionizing fronts propagate into the IGM, their expan- 
sion rate at high redshifts is limited by the mean recombina- 
tion rate in the ionized gas (rather than simply by the ioniz- 
ing luminosity). The recombination rate per unit time and vol- 
ume can be written as Q;BCHii(nHii)^, where is the recom- 
bination coefficient of neutral hydrogen to its excited states 
(= 2.6 X 10"" cm^ s"i at r = lO'^K), (nnn) is the mean num- 
ber density of ionized hydrogen, and Chii = («hii)/("hii}^ is 
the mean clumping factor of ionized gas. The clumping factor 
can be estimated by numerical simulations (Gnedin & Ostriker 
1997), and in semi-analytic models (Madau et al. 1999; Chiu 
& Ostriker 2000; Miralda-Escude et al. 2000; Haiman, Abel & 
Madau 2001; Benson et al. 2001; Wyithe & Loeb 2003). It is 
expected that the clumping factor rises from a value near unity 
to a few tens, as structures go increasingly nonlinear, between 
redshift redshifts z. = 10-20. Numerical simulations predict 
Chii = 1 — 10atz>10, but at high redshifts, where numeri- 
cal simulations cannot resolve the Jeans mass, and the effec- 
tive clumping can be increased by fluctuations on small scales 
(Haiman, Abel & Madau 2000), to Chh « 40 at z ~ 10. Finally, 
we note that whenever recombinations are fast enough to set 
the expansion rate of the ionization fronts, which is indeed the 
case in the models we consider below, it is essentially the linear 
combination c/Chh that determines the evolution of the total 
ionized volume fraction. 

2.5. Feedback Effects 

Several feedback effects are likely to be important for reion- 
ization. There can be significant internal feedback in or near 
each ionizing source, due to the presence of supernovae (Fer- 
rara 1998), or of the radiation field (Omukai & Nishi 1999; Ri- 
cotti, Gnedin & ShuU 2002), on the local H2 chemistry. The net 
sign of these effects is difficult compute, as it depends on the 
source properties and spectra, and on the detailed density distri- 
bution internal and near to the sources. For practical purposes 
of computing the global reionization history, we may, however, 
think of any internal feedback effect as regulating the efficiency 
parameter e* defined above. 

Since the universe is optically thin at soft UV (<13.6eV), and 
soft X-ray ( ^ IkeV) photon energies, radiation from the ear- 
Uest Type II halos can build up global backgrounds in these 
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bands, and provide prompt external feedback on the forma- 
tion of subsequent structures, which are easier to follow. In 
particular, soft UV photons can photodissociate H2 molecules, 
while X-rays can promote H2 formation (Haiman, Abel & Rees 
2000). The net sign of this feedback depends on the spectrum 
(with a critical X-ray/UV energy flux ratio of a few percent) 
of the earliest sources (Haiman, Abel & Rees 2000; Ciardi, 
Ferrara & Abel 2000; Machacek, Bryan & Abel 2001, 2003; 
Glover & Brandt 2003). 

A second type of important feedback is that photo-ionized 
regions are photo-heated to a temperature of ^ lO'^K, with a 
corresponding increase in the Jeans mass in these regions. As 
a result, gas from existing Type II halos is photo-evaporated 
(Shapiro et al. 1998; Barkana & Loeb 1999), and if the gas is 
kept photoionized, then Type la halos are subsequently prohib- 
ited from collecting gas at their centers. Gas infall and cooling 
in the ionized zones is suppressed in halos with virial temper- 
atures up to - 2 X 10^ K (Efstathiou 1992; Thoul & Weinberg 
1995; Navarro & Steinmetz 1997). Type la dark matter halos 
that form in already ionized regions are therefore effectively ex- 
cluded from contributing to reionization; in the ionized regions, 
new ionizing sources will have to wait until Type lb halos build 
up further down in the hierarchy of structure formation. 

3. MODELS OF REIONIZATION 

In this section, we describe our semi-analytical model of the 
reionization process, incorporating the above effects in the sim- 
plest possible way. In these models, we track the total volume 
of ionized regions, assuming that ionizing sources are located 
inside virialized dark matter halos. Each source creates an ion- 
ized Stromgren region, which expands into the IGM at a rate 
dictated by the source luminosity, and by the background IGM 
density and clumping factor. 

We adopt the mass function of virialized dark matter halos 
from the N-body simulations of Jenkins et al. (2001). The 
comoving number density {dn / dM)dM of halos at redshift z 
with mass M ± dM /2 is given by 



dn 
dM 



(z,M) = 0.315 



1 dcr, 



M 



M gm dM 



exp 



-|0.61-log(£),aM; 



,3.8 



(1) 

where cjm is the r.m.s. density fluctuation, computed on mass- 
scale M from the present-day Unear power spectrum (Eisen- 
stein & Hu 1999), is the Unear growth function, and po is 
the present-day total mass density. 

The gas inside each viriaUzed object is shock-heated to the 
virial temperature, 

/ M /o„\ 1/3 / u \2/3 
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where /i is the mean molecular weight (p, = 0.60 for ionized gas 
collecting in Type Ia,b halos, and /u = 1 .22 for neutral gas of pri- 
mordial composition, relevant for gas collecting in the shallow 
potentials of Type II halos). 

In order to address the two global radiative feedback effects 
described in § 2.5 above, we separately define, and track, the 
total fraction of the mass in the universe that is condensed into 
Type II, la, and lb halos: 
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Fig. 1 . — The evolution of the collapsed fraction of all baryons in halos in 
three different virial temperature ranges. The three different virial temperature 
ranges represent halos that can form ionizing sources (1) only in the neutral 
regions of the IGM, and only if H2 molecules are present, (2) only in the 
neutral IGM regions, irrespective of the H2 abundance, and (3), irrespective of 
the H2 abundance, and both in the ionized or neutral phases of the IGM. 



and 



^colflb' 



{z) = - r dM(^\M, (5) 

where the halo masses Mu = Mu(7ii,z), Mi^ = Mi^(Tia,z), and 
Mib = A/ib(rib,z) are obtained from equation (2) for the temper- 
atures Til = lOOK, Tia = 10"* K and Tib = 2 X 10^ K. Our results 
below are insensitive to the upper limit in equation (5). The 
collapsed fractions are shown in Figure 1. The time deriva- 
tives of equations (3)-(5), in turn, represent the rates at which 
mass is condensing into nonlinear halos. At low redshifts, the 
derivatives can turn negative; which represents the decrease in 
the number of halos as they merge into larger systems. This 
does not significantly effect our results at the high redshifts we 
consider z ^ 6. 

We assume that each halo represents a single cosmological 
"reionizing source", which drives an expanding ionized region 
into the IGM. For each initial source redshift, we solve the 
equation of motion for the ionization front, /?,, 



dt ' 47r(n/f) 



-C}m{nH)aBR], 



(6) 



where Hiz) is the Hubble constant at z, q;b = 2.6 x 10~'^ cm^ s"' 
is the recombination coefficient of neutral hydrogen to its ex- 
cited states at T = lO'^K, and Cnn is the mean clumping factor of 
ionized gas within as discussed in § 2.4 above. The first term 
in equation (6) accounts for the Hubble expansion, the second 
term accounts for the ionizations by newly produced photons, 
and the last term describes recombinations (Shapiro & Giroux 
1987; Haiman & Loeb 1997). 

In equation (6), My represents the emission rate of ionizing 
photons from a single halo. For a stellar population with a 
Salpeter IMF (initial mass function) undergoing a burst of star 
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formation at a metallicity equal to 2% of the solar value (Lei- 
therer et al. 1999), we find that the rate is well approximated 

by 

N(t)-!^° (f<10^-^yr) 



(t > 10^-5 yr), 



(7) 



where A^o = 3.7 x 10 s"' Mq (per Mq of stellar mass). Over 
the lifetime of the population, this produces w 4000 ionizing 
photons per stellar proton. In practice, the injection of ioniz- 
ing photons is nearly instantaneous compared to the expansion 
timescale even at redshift z ^ 20, and our results are insensitive 
to the adopted lightcurve (an assumption that remains justified 
for an IMF biased to more massive stars). 

In order to investigate the various physical effects outlined 
in § 2, we define separate star-formation efficiencies, eu, eia, 
and eib for Type II, I, and lb halos, respectively. In order to 
allow for a negative radiative feedback by the soft UV back- 
ground on the H2 abundance, we consider a redshift Zm below 
which en = 0, so that Type II halos are excluded from contribut- 
ing to reionization. In reality, rather than dropping to zero, the 
efficiency in Type II halos should roughly remain at a constant 
"steady state" value, set by the level of starformation that can 
just be maintained against the H2 photo-dissociation implied 
by this starformation (Haiman, Abel & Rees 2000). Retaining 
this level of starformation in our models would have a relatively 
small effect on our results. 

For the sake of simplicity, we adopt a constant clumping fac- 
tor Chii. While the clumping factor does evolve with redshift, at 
z ^ 10, this could be absorbed as a factor of ^ 10 decrease in the 
efficiencies e down to this redshift. Including the evolution of 
the clumping factor with redshift would be necessary in order to 
model the evolution of the neutral fraction at 5.5 ^ z ^ 6.5, nec- 
essary to interpret its sharp rise in the observed high-redshift 
quasar spectra. However, semi-analytical models (Haiman & 
Loeb 1997, 1998) and numerical simulations (Gnedin & Os- 
triker 1997; Gnedin 2001; Razoumov et al. 2002; ) have shown 
that even though the entire process of reionization can last for 
an extended period, the rapid evolution in the neutral fraction, 
caused by the sudden buildup of the radiation background once 
the ionized regions percolate, follows the percolation epoch 
within a short redshift interval Az < 1 . We will require only 
that our model indicate a "second reionization" epoch at z 7, 
and defer detailed modeling of the neutral fraction to a future 
paper 

Combining the abundance of the ionizing sources (eq. 1) with 
the ionized volume associated with each source (eq. 6), we can 
compute the volume filling fraction of the ionized regions as a 
function of redshift. The ionized volume associated with a halo 
of mass M scales linearly with its mass and with the efficiency 
factor e), namely yHii(Zon,z,M) = {AiiRl/3) = eM x Vuii(Zon,z), 
where VHu{Zon,z) is the ionized volume per unit mass and unit 
efficiency, and Ri is the solution of equation (6). The HII filling 
factor pHiiiz) is then obtained from the collapsed gas fractions 
as 

Fun (z) = p,{z) rjz'(eib^%ii^(z')+[l-ftii(z')] x 



dz 

eia ; (z ) + eii ; (z ) 



dz 



dz 



Vmi(z,z), 



(8) 



where pb = f^bPcrft is the average baryonic density and pcrit is 
the critical density of the universe. In the second term on the 
right hand side of equation (8), we have explicitly included a 




Fig. 2. — The evolution of the ionized fraction of all baryons. The thick 
solid curve corresponds to our fiducial model described in § 4.1. The other 
curves indicate variations in the importance of H2 photodissociation feedback, 
from no feedback (top dot-dashed curve), to strong feedback that effectively 
excludes Type II halos from contributing to reionization (bottom dot-dashed 
curve). The pair of dotted curves bracket a range At = ±0.01 around the 
fiducial model; the short and long-dashed curves correspond to reductions in 
T by At = 0.03 and 0.06, respectively (see § 4.2 for discussion). 



factor (1 -Full), which takes into account the fact that new ion- 
izing sources should appear in Type II and Type la halos only 
in regions which have not yet been ionized. 

Equations (1-8) determine the evolution of the ionized frac- 
tion. Note that equation (8) is implicit, and, in practice, has 
to be solved numerically by starting at high redshift (z = 60), 
and taking sufficiently small negative redshift steps (we found 
Az ~ 0.001 to be adequate in most of the models we studied). 
In summary, our model has five parameters: the clumping fac- 
tor Chii, the overall efficiencies, en, eia, ejb, and the redshift 
Zuv at which H2 dissociative feedback sets in. Since variations 
in the clumping factor are nearly degenerate with changes in 
the efficiency parameters, in all of our models we fix Chii =10 
(alternatively, adopting a redshift-dependent clumping factor 
from the works discussed in § 2.4 would only case a "tilt" in 
the reionization histories we derive below). We therefore have 
four "free" parameters. Current CMB constraints are not suffi- 
ciently accurate to converge on a "best-fit" model, and our main 
purpose here is to be able to study the various physical effects 
outlined in § 2 above. Nevertheless, the hope (as we attempt to 
demonstrate in § 5 below) is that eventually, CMB polarization 
data will provide constraints on parameters of models similar 
to the one described here, while simultaneously delivering con- 
straints on cosmological parameters. 

4. POSSIBLE REIONIZATION HISTORIES 

In this section, we present a suite of models, designed to in- 
vestigate the various effects summarized in § 2 above. Note 
that more realistic calculations would explicitly have to couple 
the star formation efficiency with radiative feedback processes 
(as attempted in part by Haiman, Abel & Rees 2000, Glover & 
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Fig. 3. — The evolution of the ionized fraction of all baryons. with the thick 
solid curve showing the fiducial model as in Figure 2. The other curves indicate 
variations in the efficiency parameter en in Type II halos. As in Figure 2, 
the pairs of dotted, short-dashed, and long-dashed curves bracket a range of 
Ar = ±0.01, 0.03, and 0.06 ai'ound the fiducial model. The required changes 
in the efficiencies are described in the text. The bottom dot-dashed curve 
assumes Type II halos do not form stars, and the top dot-dashed curve assumes 
that there is no H2 photodissociation feedback (see § 4.3 for discussion). 



Brandt 2003, and Cen 2003). In the present models, we explic- 
itly track only the ionized fraction of hydrogen. For the purpose 
of computing the electron scattering optical depth, we assume 
that the HII and Hell fractions are identical, so that the extra 
electrons from Hell increase the opacity by a factor of 1.08. A 
Hell to Helll transition appears to be taking place at the low 
redshift z - 3 (Miralda-Escude 2001). While some Helll could 
be present at high redshifts, for realistic spectral slopes of the 
ionizing sources, the extra electrons from Helll regions add cor- 
rections to T at the <^ 10% level. 

4.1. Fiducial Model 

As discussed above, in order to match inferences from the 
z < 6 quasar spectra, we require that a percolation takes place 
at redshift z ~ 7. A second motivation for requiring a percola- 
tion near z ^ 7 is the high observed temperature of the IGM at 
z ~ 4: if the percolation occurred too early, the IGM would 
cool to temperatures that may be too low by z ~ 4 (Hui & 
Haiman 2003). We find that this effectively fixes the efficien- 
cies eia = ejb = 80 (if they are postulated to have the same value). 
Accordingly, in our fiducial model, we adopt these values. Note 
that these values are quite reasonable: for a normal stellar pop- 
ulation, N~f = 4000, and we may break down the efficiency into 
/* =0.2 and /esc = 0.1. With eia and eib = 80 fixed, we have 
only Zuv and en to vary. In our fiducial model. Type II halos 
host metal-free stars from a small fraction (/, = 0.005) of the 
available gas, but the stars are, on average, 10 times more effi- 
cient ionizing photon producers (A^.^ = 40000), and all of their 
ionizing radiation escapes into the IGM (/esc = 1). This results 
in the overall efficiency of eu = 200. We then finally adopt a 
redshift Zuv = 17 for the onset of H2-photodestruction, which 



Fig. 4. — This figure assumes that H2 feedback effectively shuts off starfor- 
mation in all Type 11 halos, and therefore Type la halos form metal free stars, 
but Type lb halos form normal stars. The thick solid curve shows a new fidu- 
cial model that has the r = 0.17. This requires increasing the efficiencies in 
Type la halos by a factor of 60 relative to the fiducial models in Figures 2 and 
3. The pairs of dotted, short-dashed, and long-dashed curves bracket a range 
of At = ±0.01, 0.03, and 0.06 around this new fiducial model. Note we keep 
the efficiency in Type II halos fixed, and percolation in all models is achieved 
only at z ~ 6. This can lead to consistency with inferences from z ~ 6 SDSS 
quasar spectra (see § 4.4 for discussion). 



allows Type II halos to ionize ^ half of the volume of the IGM. 
By construction, this model produces an optical depth r = 0.17, 
the central value in the range allowed by WMAP. In what fol- 
lows, we separately vary each parameter - this will explicitly 
clarify the dependence of the reionization history and of r on 
each of our parameters. 

4.2. When does H2 dissociation feedback set in? 

The redshift at which any H2 dissociation sets in is deter- 
mined by the spectra of the light sources in Type II halos, and 
the intensity of the X-ray and soft UV backgrounds. According 
to Haiman, Abel & Rees (2000), the choice of Zuv = 17, result- 
ing in Type II halos ionizing ^half of the IGM volume, would 
be consistent with an X-ray/UV background energy flux ratio 
of a few percent. We here simply parameterize this feedback 
with the redshift of its onset, Zuv The reionization history in 
our fiducial model is shown by the thick solid curve in Figure 2. 
The other curves in this figure show variations in Zuv away from 
the fiducial model, which correspond to changes Ar = 0.01, 
0.03, and 0.06. The level of the constraint that will be provided 
by Planck corresponds approximately to Ar = 0.01 (although 
Planck will have information beyond a measurement of r, see 
discussion in § 5 below and in paper II). As the dotted curves in 
the figure show, at the level of the Planck constraint, Zuv could 
be determined to a precision of Azuv ~ ±1.5. Note that even 
in the absence of any H2 destruction (upper dot-dashed curve), 
the opacity can only be increased to r 0.185. 

4.3. How different is the overall efficiency in Type II halos? 
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In our fiducial model, we have adopted a ratio cn/eia = 2.5, 
determined by the transition from metal-free to normal stars. In 
Figure 3, we consider variations in eu. We show results for val- 
ues of en that result in changes of ±At = 0.01, 0.03, and 0.06 
(dotted, short-dashed, and long-dashed curves, respectively). 
In each case, we adjust the value of Zuv so that the Type II halos 
always ionize ^ half of the volume (this approximately mim- 
ics fixing the radiation intensity at which H2 feedback sets in). 
The value of en for these six curves, from bottom to top, are 
en = 20,80,120,280,520, and 1400. The last value still rep- 
resents a relatively modest factor of ^ 7 increase in the com- 
bined efficiencies, and is within the range of uncertainties. The 
Planck-level constraints represent ^ 50% changes in en. The 
bottom dot-dashed curve shows the case of no starformation in 
Type II halos (en = 0; this is equivalent to Zuv ^ 35, shown by 
the bottom curve in Figure 2). The top dot-dashed curve shows 
the highest efficiency case, but assuming no H2 dissociation. In 
this model, r = 0.29, which would be ruled out at the ~ 3(t level 
by the current WMAP results. 

4.4. Are H2 molecules required to produce large r ? 

In Figure 4, we assume that H2 feedback effectively shuts 
off starformation in all Type II halos. As a result. Type la ha- 
los form metal-free stars, and, motivated by the arguments in 
§ 2.3, we assume that the transition to normal stellar popula- 
tions took place in Type lb halos (and we set eib = 80). We then 
revise our fiducial model and find that eia = 4800 is required to 
produce t = 0.17, implying an overall efficiency difference of a 
factor of 60 between Type la and lb halos (shown by the solid 
curve). The main conclusion from Figure 4 is that r = 0.17 
can be achieved without H2 molecules, but this requires a boost 
in efficiencies by this factor This would require, in turn, in 
addition to the maximum factor of ^ 20 caused by switching 
from metal-free to normal stars, that the escape fraction or star- 
formation efficiency in Type la halos was a factor of 3 higher 
(or else that the clumping factor decreases by the same fac- 
tor from 15 to 10). Another important feature apparent 
from Figure 4 is that the reionization history is qualitatively dif- 
ferent from those shown in Figures 2 and 3. As we discuss in 
§ 5 below, future CMB polarization measurements can distin- 
guish between the two scenarios at high significance. As in the 
previous figures, we also show the ±At = 0.01, 0.03, and 0.06 
"contours", which correspond (bottom to top curves) to eia = 
480, 1600, 3200, 8000, 16000,48000. Planck-level constraints 
on eja = are at the ~ 50% level. It is worth pointing out that in 
the models shown in Figure 4, Type la sources are self-limiting 
by photo-heating feedback (shown by the flattening of the ion- 
ization histories below z <^ 12) and it is the Type lb sources that 
eventually cause percolation at z 6.5 in all models, regardless 
of the efficiency in Type la sources. 

4.5. Constraints on efficiencies from the percolation epoch 

In our models. Type lb halos form last, and must be respon- 
sible for the eventual percolation of the ionized regions. As ar- 
gued above, to be consistent with inferences from z ^ 6 quasar 
spectra, this percolation epoch cannot be at z ^ 7. It is interest- 
ing to ask, then: what is the allowed efficiency range in Type lb 
halos? To address this question, we consider the fiducial model 
defined in § 4. 1 , and we vary eja = eib . The resulting reionization 
histories for variations that cause ±Ar = 0.01, 0.03, and 0.06 
are shown in Figure 5. The figure shows that in high-efficiency 
models (corresponding to At ^ 0.03, or e/i ^ 1200), percola- 




FlG. 5. — The evolution of the ionized fraction of all baryons. with the 
thick solid curve showing the fiducial model as in Figure 2. The other curves 
indicate variations in the star formation efficiency in Type la and lb halos . 
The pairs of dotted, short-dashed, and long-dashed curves bracket a range of 
At it 0.01, 0.03, and 0.06 around the fiducial model, as in Figure 3 (see § 4.5 
for discussion). 

tion occurs at too high redshifts. Similarly, models with too low 
efficiencies (At <, -0.01, or ejb ^ 30) are ruled out as they fail 
to reionize the universe by redshift z ^ 6. The low-redshift ob- 
servations by themselves therefore constrain the efficiencies in 
Type lb halos to within a range of a factor of ^ 30. 

4.6. The running spectral index WMAP model 

When combined with other large scale structure data, the 
WMAP results favor, at the ~ 2a level, a non-scale-invariant 
primordial power spectrum with a running spectral index (Spergel 
et al. 2003). This model has much less small scale power than 
the best-fit WMAP-only model. It is interesting to ask whether 
T = 0. 17 can still be achieved in the running index model with- 
out H2 molecular cooling (Spergel et al. 2003). To address this 
question, we consider the fiducial model defined in § 4.4, but 
use the cosmological parameters listed in Table 10 of Spergel 
et al. (2003). We then increase eja until r = 0.17 is achieved. 
In Figure 6, we show the reionization history in this model, 
which requires boosting the efficiencies by a factor of 3000, to 
eia = 2.5 X 10^. This model must unrealistically assume that 
all the gas turns into stars (and gets recycled ^3 times), all of 
the ionizing radiation escapes into the IGM, the clumping fac- 
tor is Cen = 1, and adopt the maximum factor of ^ 20 in the 
increase in ionizing photon production rates in massive, metal- 
free stars. The dotted, short-dashed, and long-dashed curves 
assume eja = 80000, 20000, and 9000, corresponding to decreas- 
ing T by At = -0.01,-0.03, and -0.06, respectively, and show- 
ing that the lack of molecules could still be consistent with the 
lower range of r allowed in the running index model. Allowing 
for Type II halos (assuming no H2 feedback) we find that a rel- 
atively more realistic e// = 4000 is required to produce t = 0. 17. 

4.7. Summary 
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Fig. 6. — This figure replaces the best fit power law WMAP cosmological 
model by the best fit model that allows a running spectral index and combines 
WMAP and other datasets (Table 10 in Spergel et al. 2003). In this model, 
small-scale power is significantly reduced. The thick solid curve shows a 
new fiducial model that assumes no star-formation in Type 11 halos, and has 
T = 0.17. This requires boosting the efficiency by an unrealistic factor of 3000 
relative to our fiducial model in the power-law cosmology. The dotted, short- 
dashed, and long-dashed curves correspond to lower efficiencies, and reduced 
values of t by At = —0.01, —0.03, and -0.06, respectively (see § 4.6 for dis- 
cussion). 



The physics of reionization is rich in features that can natu- 
rally lead to distinctive ionization histories. These features can 
arise because of (1) the different types of coolants in Type II 
and Type I halos, (2) the different response of different halos 
to radiative feedback on the H2 chemistry, and to photoioniza- 
tion feedback on gas infall, and (3) the different efficiencies ex- 
pected in metal-free and normal stellar populations. We generi- 
cally find that no simple model fits both the WMAP data and the 
constraint from z ^ 6 quasar spectra requiring that an episode 
of percolation occur near z < 7. Nevertheless, a range of inter- 
esting, physically motivated reionization histories is allowed by 
current data. We find that the ionized fraction Fuii(z) in these 
models evolves monotonically with redshift, except under the 
assumption of a uniformly increasing metal-enrichment of the 
IGM, and a sudden transition from metal-free to "normal" stars 
at some redshift Za- (see discussion in next section, and also 
Wyithe & Loeb 2003 and Cen 2003). 

5. CONSTRAINTS FROM FUTURE CMB ANISOTROPY 
MEASUREMENTS 

All models shown in Figures 2-6 in the previous section are 
consistent with the optical depth t = 0.17 ±0.04 measured in 
the current WMAP data. These figures also give an indication 
of constraints that will be available in the future on reioniza- 
tion models - by showing variations in each model parameter 
that cause changes as small as At = ±0.01 in the electron scat- 
tering optical depth. While the constraints on single param- 
eters are impressive at this At, considering one parameter at 
a time will clearly be only the first crude step in attempts to 
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Fig. 7. — This figure addresses the hypothesis that a transition from metal 
free to normal star formation takes place abruptly at a fixed redshift 

Zcr, accom- 
panied by a sudden decrease by a factor of in the ionizing photon production 
rates. The thick solid curve shows a model that has t = 0.17; this is achieved 
by Zci = 14 and fa = 20. The pairs of dotted, short-dashed, and long-dashed 
curves bracket a range of At it 0.01, 0.03, and 0.06 around this model by 
changing Za- 



understand the reionization history. Eventually, a full degen- 
eracy study, involving simultaneous variations of all reioniza- 
tion model parameters, together with cosmological parameters, 
will be needed. While such a study is not yet called for by the 
WMAP data, our results indicate that this will be necessary for 
Planck (see paper II for in-depth discussion). 

Such future studies will utilize information in the CMB po- 
larization and temperature power spectra beyond the value of 
the optical depth r. Indeed, the lack of our a-priori knowledge 
of the reionization history would bias estimates of t itself, and 
would be the limiting factor in determining r from Planck data 
(see paper II). Here we consider the question: how useful is 
the information beyond the measurement of r in constraining 
reionization histories? 

To address this question, we here consider three different 
models that have the same value of t = 0.17. Two of these 
models are the fiducial models defined in § 4. 1 and § 4.4, and 
shown by the solid curves in Figures 2 and 4; hereafter referred 
to as models (A) and (B). The electron scattering in model (A) 
extends to redshifts higher than in model (B), and, as a result, 
will produce polarization anisotropy on smaller angular scales 
(larger spherical harmonic index £). For sake of illustration, we 
consider a third model, similar to those proposed by Wyithe & 
Loeb (2003) and Cen (2003). This is a variant of model (A), 
in which we assume that a transition from metal free to nor- 
mal stellar populations takes place abruptly at a fixed redshift 
Zcr, accompanied by a sudden decrease by a factor of f„ in the 
ionizing photon production rates. In terms of our model param- 
eters, £// = e/a = eib = 80 at z < Zci- and e// = e/„ = e/i = 80/cr at 
z > Zcr- We find that Za- = 14 and f„ = 20 results in t = 0.17; 
this model is shown in Figure 7 (together with models with 
Zcr = 10,11,13,15,17, and 21, which bracket the Ar±0.01, 
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Fig. 8. — The large-angle E-mode polarization anisotropies in three dif- 
ferent models. The dashed, dotted, and solid curves correspond to the solid 
curves in Figures 2, 4, and 7, respectively. All three models produce the same 
electron scattering opacity t = 0.17. Cosmic variance error bars are shown 
around the solid curve, in bins of A£ = 1 . The three curves are statistically 
distinguishable in a cosmic variance limited experiment at > 4(t significance. 
Planck can only distinguish the dashed and the solid curves at ~ 1.6(t, but can 
distinguish the other two pairs of models > 3cr (see Table 2). 
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Fig. 9. — The same three models as in Figure 8, but the temperature- 
polarization cross-correlation is shown. The three models are nearly indis- 
tinguishable. 



0.03, 0.06 range around r = 0.17). Our motivation to reconsider 
this model, hereafter referred to as model (C), is the distinctly 
different shape of the reionization history. 



We modified CMBFast^ (Seljak & Zaldarriaga 1996) to use 
the ionization histories from models (A), (B) and (C) to gener- 
ate temperature, polarization, and cross-correlation anisotropy 
power spectra, CTe,CEe and Cce respectively. The resulting EE 
and TE power spectra are shown by the dashed, dotted, and 
solid curves in Figures 8 and 9, respectively. In Figure 8, we 
show cosmic variance error bars in bins of A£ = 1 around model 
(C), and in Figure 9, we overlay the WMAP data. The models 
are nearly indistinguishable in the TE power spectra, but are 
clearly different in the EE spectrum: model (A) is biased to 
lower £ than model (B), and model (C) has excess power at 
small angular scales. 

We next calculate the likelihood between pairs of models, us- 
ing all three (EE, TE, and TT) power spectra, and a Ax^ statis- 
tic defined in paper 11 (see also Kaplinghat et al. 2003). As 
apparent from Figures 8 and 9, nearly all of the distinguishing 
power comes from the E-mode polarization anisotropy alone. 
In Table 2, we show the significance at which pairs of mod- 
els can be distinguished in the future by (1) 2-year of WMAP 
data (2) Planck, and (3) an "ultimate" cosmic variance limited 
experiment. The details of folding the instrumental characteris- 
tics of WMAP and Planck into our calculations are described in 
paper 11. 

The three curves are statistically distinguishable in a cosmic 
variance limited experiment at > 4a significance. Planck can 
only distinguish models (A) and (C) at 1 .6cr, but can distin- 
guish the other two pairs of models at > 3a significance (see 
Table 2). It is also interesting to ask how poor a fit a "one-step" 
reionization model would be to real data realized from a more 
complex history? This issue is discussed in detail in paper 11; 
for sake of completeness, we quote here only the difference of 
model (A) from one-step reionization models. We find likeli- 
hoods of Ax^ of 0.3, 16, and 60, for 2-yr WMAP data, Planck, 
and a cosmic variance limited experiment, respectively. 

In summary, our results suggest that for the physically moti- 
vated reionization histories discussed in this paper, the Planck 
satellite will be able to rule out both a simple 1-step reion- 
ization model, and also several "wrong" physically motivated 
models, at high statistical significance. 

6. CAN JWST DETECT THE SOURCES OF REIONIZATION? 

A broad conclusion that can be drawn from the WMAP re- 
sults, reinforced in this paper, is that UV sources exist at very 
high redshifts (z ^ 20). This, by itself, is a significant new dis- 
covery, and it raises the question of whether we will be able 
to directly detect these sources. Here we consider the capabil- 
ities of the James Webb Space Telescope (JWST): what is the 
"limiting redshift" out to which JWST will be able to detect the 
sources of reionization? We here answer this question with a 
simple calculation, as follows. Consider a halo of mass M at 
redshift z- The number of halos A^(> M) with mass above M at 
redshift z±dz/2 that could be visible within a a soUd angle Ai7 
and redshift range Az is given by 



A^(> M) = A17 X Az X /o, 



dV 
dzdfl 



(z) / dM' 



dn 
dM' 



(9) 



where dV /dzdH. is the cosmological volume element, dn/dM 
is the halo mass function from Jenkins et al. (2001), and /on 
is the fraction of the time the halo spends at a luminosity de- 
tectable by JWST. We assume the sources live for f, = 3 x 10^ 
years (approximately the main sequence lifetime for massive 

^ available at http://www.physics.nyu.edu/matiasz/CMBFAST/onhne.html 
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stars, both normal and metal-free, and a conservative estimate 
for the bright phase of mini-quasar activity). In this case, /on « 
tslt{z) « 0.006-0.02 at redshifts z = 10-20 (where t{z) is the 
age of the universe at redshift z). 

The sensitivity expected to be reached in a 5 x 10"*second ex- 
posure "yWSr Deep Field", with the current 6m diameter tele- 
scope design, in the 0.6-5/um band (S/N=5), with the Near 
Infrared Camera (NIRCam) is ~ InJy. We require that there 
should be at least 100 sources visible down to this depth in the 
AJl = 1 arcmin- field of view of NIRCam, and in Az ?» 2. To 
predict the luminosities, consistent with our reionization mod- 
els, we make the assumption that in a halo of mass M, ~ 10% 
of the gas turns into stars with a normal Salpeter IMF. We use 
the population synthesis models of Leitherer et al. (1999) to 
predict the observed spectra (at a metallicity of Z = 0.02). We 
ignore wavelengths shorter than 1215(1 -i-z)A, because of the 
strong GP trough expected at high redshifts, and compute the 
average flux in the unobscured fraction of the 0.6 — 5/um wave- 
length range range. 

Under these assumptions, we find that a halo of mass 2 x 
10^ Mq at redshift z = 14 would have an average flux of ^ InJy 
in the unobscured 1.8-5/im range, detectable with NIRCam. 
Such a halo represents a 3(t object at z = 14 in the WMAP cos- 
mology, and with our derived duty-cycle of /on = 0.01 there 
should be 100 such sources detectable at any given time in 
10 arcmin^. This halo would have a virial temperature of ^ 
2 X lO'* K, justifying the assumed star formation efficiency of 
10% (larger than is expected in Type II halos). As we argued 
above in § 2.3, it is possible that most stars are still metal-free 
in Type la halos at these high redshifts. Although metal free 
stars produce more ionizing photons than normal stars, owing 
to their high effective temperatures, they are fainter at the JWST 
wavelengths than normal stars. This, however, is likely offset 
by the steep IMF expected from a metal-free stellar popula- 
tion. We conclude that JWST could directly detect the sources 
of reionization, out to redshifts as high as z ^ 14; these sources 
would be sufficiently bright to allow follow-up spectroscopy 
with NIRSpec on JWST. The sources could also be discovered 
in line emission, either in the Hf3 line (Oh 1999), or, in case 
of metal-free stellar populations, in the He 1640A, and 4686 A 
lines out to redshifts z ^ 20 with the mid-infrared instrument 
(Oh, Haiman & Rees 2001; Tumlinson, Giroux & ShuU 2001). 

7. CONCLUSIONS 

The determination by WMAP of a high optical depth to Thom- 
son scattering implies that significant reionization took place at 
redshifts z ^ 20, and has provided our first glimpse into the 
earliest structures that exist at these distant epochs. In this pa- 
per, we used physically motivated semi-analytical reionization 
models to quantify the implications of the WMAP result. The 
high value of r, by itself, requires that the earliest light sources 
have a surprisingly high efficiency e of injecting ionizing pho- 
tons into the IGM. In the running spectral index cosmological 
model favored by the combination of WMAP with other large 
scale structure data (Spergel et al. 2003), H2 molecules are 
required to form efficiently in halos with virial temperatures 
Tvir < lO'^K at high redshifts to produce r = 0.17. 

Our most interesting conclusions arise when we require that 
our models produce (1) a high optical depth r = 0.17, and (2) 
a percolation epoch at z ^ 6-7. The latter is required by most 
interpretations of the rapid evolution in the neutral fraction in 
the IGM at z = 5 - 6 (Songaila & Cowie 2002; Becker et al. 



2001), and possibly by the high temperature of the IGM in- 
ferred from the Lya forest at z ~ 4 (Hui & Haiman 2003). We 
genericaUy find that no simple model can simultaneously sat- 
isfy both of the above requirements. Our results indicate that 
either H2 molecules had to form efficiently at z ^ 20 and al- 
low UV sources in halos with virial temperatures Jvir < 10'*K 
to contribute substantially to reionization, or else the ionizing 
photon production efficiency in halos with T^ij > IQ^K had to 
decrease by a factor of ^ 30 between z ~ 20 and z ~ 6. The 
first of these two options would have important implications 
on radiative feedback processes that determine the global H2 
chemistry. The second option would naturally be explained by 
the presence of metal-free stars in the earliest halos, which are 
then replaced by normal stellar populations in larger halos that 
collapse at later epochs. 

The physics of reionization is rich in features that can nat- 
urally lead to distinctive ionization histories. These features 
can arise because of (1) the different types of coolants in halos 
with virial temperatures above and below ~ lO'^K, (2) the dif- 
ferent response of different halos to radiative feedback on the 
H2 chemistry, and to photoionization feedback on gas infall, 
and (3) the different properties of metal-free and normal stel- 
lar populations. A range of physically motivated reionization 
histories is allowed by current data. We find that the ionization 
history in our models generally evolves monotonicaUy with red- 
shift, in difference from "double reionization" that results if the 
metallicity was increasing uniformly in the IGM, triggering a 
sharp transition at a fixed redshift from metal-free to "normal" 
stars. Nevertheless, we find that the evolution of the ionized 
fractions have sufficiently conspicuous features for Planck to 
provide tight statistical constraints on reionization model pa- 
rameters, and to elucidate much of the physics at the end of 
the Dark Ages. The sources responsible for the high optical 
depth discovered by WMAP should be directly detectable out to 
z ~ 15 by the James Webb Space Telescope. 

The WMAP result has opened a new window into studies of 
the first structures at the end of the cosmological dark ages, 
and has raised the hopes for further interesting findings from 
all potential probes of the high-redshift universe; including fu- 
ture CMB polarization anisotropy studies, spectroscopic obser- 
vations of quasars and Lya emitting galaxies at z > 6, radio 
probes of the redshifted 21cm line of neutral hydrogen, as well 
as direct detections of the sources of reionization and their end- 
products as supemovae and gamma ray bursts. 

We thank Manoj Kaplinghat, Lloyd Knox, and David Spergel 
for stimulating discussions. GPH is supported by the WM. 
Keck foundation. We thank Uros Seljak and Matias Zaldarriaga 
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Table 1 

Collapse redshifts and masses of 1, 2 and 3o- dark matter halos with virial temperatures of Tvir = 100 and 10'*K, 

REFLECTING THE MINIMUM GAS TEMPERATURES REQUIRED FOR COOLING BY H2 AND NEUTRAL ATOMIC H, RESPECTIVELY. 



rvir(K) 


V 


Zcoll 


Mhalo(M0) 


100 


1 


7 


2x 10^ 




2 


16 


5x lO'' 




3 


26 


3 X lO'' 


lO'' 


1 


3.5 


4x 10** 




2 


9 


1 X 10* 




3 


15 


6x 10' 



Table 2 

Constraints (Ax^) from future CMB polarization experiments on three models with the same t = 0.17. 



Model pair 


2-yr WMAP 


Planck 


Cosmic variance 


AB 


0.2 


9 


25 


AC 


0.1 


2.6 


18 


EC 


0.3 


13 


34 


A vs one-step 


0.3 


16 


60 



